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The concept of Space Manifold Dynamics is a new method of space research. We have applied it along with the basic 
idea of the method of Ott, Grebogi and York (OGY method) to stabilize the motion of a spacecraft around the triangular 
Lagrange point L 5 of the Earth-Moon system. We have determined the escape rate of the trajectories in the general three- 
and four-body problem and estimated the average lifetime of the particles. Integrating the two models we mapped in detail 
the phase space around theLs point of the Earth-Moon system. Using the phase space portrait our next goal was to apply 
a modified OGY method to keep a spacecraft close to the vicinity of L 5 . We modified the equation of motions with the 
addition of a time dependent force to the motion of the spacecraft. In our orbit-keeping procedure there are three free 
parameters: (i) the magnitude of the thrust, (fi) the start time and {Hi) the length of the control. Based on our numerical 
experiments we were able to determine possible values for these parameters and successfully apply a control phase to a 
spacecraft to keep it on orbit around L 5 . 
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1 Introduction 

The concept of Space Manifold Dynamics (SMD) encom¬ 
passes the various applications of methods to analyze and 
design spacecraft missions, ranging from the exploitation 
of libration orbits around the collinear Lagrangian points 
(denoted by Li,L 2 ,L^) to the design of optimal station¬ 
keeping and eclipse avoidance manoeuvres or the deter¬ 
mination of low energy lunar and interplanetary transfers 
(IPerozzi & Ferraz-Mello 201 Oi l. 

The classical calculations of interplanetary trajectories 
(e.g. Apollo lunar missions, the Voyagers grand tour, or the 
New Horizons mission) are based patched-conics transfers 
and using the ’’gravity assist” principle exploiting the grav¬ 
itational force of planets. The demand for new and unusual 
kinds of orbits to meet their scientific goals made it nec¬ 
essary to look for a new approach which uses the dynam¬ 
ics of manifolds. The available tools allows researchers to 
study and analyse the natural dynamics of the problem and 
provides a way to, for instance, design new ’’low-energy” 
orbits, cheap station keeping procedures and eclipse avoid¬ 
ance strategies. These tools are based on the examination of 
the stable and unstable manifolds of the system. 

The International Sun/Earth Explorer 3 (ISEE-3, 
launched on August 12, 1978) was the first artificial object 
placed at Li, proving that such a suspension between grav- 
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itational fields was possible. The Solar and Heliospheric 
Observatory (SOHO, launched on December 2, 1995) was 
the successor of the ISEE-3 ( jPunham & Earquhar 2003| . 
The SMD approach was first applied to the SOHO pro¬ 
gram (IGomez et al. 200 it . SOHO was in Li and explored 
the Sun-Earth interactions, first examining the phenomenon 
known as ’’space weather” now. After the success of these 
missions the scientific interest increased towards the bene¬ 
fits of Lagrangian points and the new method was used to 
design more and more ambitious missions that were later re¬ 
alized in the MAP, Genesis, and Herschel-Plank missions. 
Genesis was the first spacecraft whose orbit has been en¬ 
tirely designed with the new method (IHowell et al. 1998] l. 

The Lagrangian points have obvious benefits for space 
missions. Eor example, the Li point of the Earth-Sun sys¬ 
tem affords an uninterrupted view of the Sun and is cur¬ 
rently home to the SOHO. The L 2 point of the Earth-Sun 
system was the home to the WMAP spacecraft, current 
home of Planck, and future home of the James Webb Space 
Telescope. L 2 is ideal for astronomy because a spacecraft is 
close enough to readily communicate with Earth, can keep 
Sun, Earth and Moon behind the spacecraft for solar power 
and (with appropriate shielding) provides a clear view of 
deep space for our telescopes. The Li and L 2 points are 
unstable on a time scale of approximately 23 days, which 
requires satellites orbiting these positions to undergo regu¬ 
lar course and attitude corrections. It is unlikely to find any 
use for the L 3 point since it remains hidden behind the Sun 
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at all times. The idea of a hidden ”Planet-X” at the L 3 has 
been a popular topic in science fiction writing. The instabil¬ 
ity of Planet X’s orbit (on a time scale of 150 years) didn’t 
stop Hollywood from shooting movies like The Man from 
Planet X. 

The L 4 and L 5 points are home to stable orbits so long 
as the mass ratio between the two large masses exceeds 
24.96 (in our case between the Earth and the Moon it is 
81) ( Murray & Dermott 1999| l. This condition is satisfied 
for both the Earth-Sun and Earth-Moon systems, and for 
many other pairs of bodies in the solar system. In 2010 
NASA’s WISE telescope hnally confirmed the first Trojan 
asteroid (2010 TK7) around Earth’s leading Lagrange point. 

The stability property of L 4 and L 5 offers an excel¬ 
lent opportunity for spacecrafts to maintain their orbits with 
minimal fuel consumption. These points could be excellent 
locations to place space telescopes for astronomical obser¬ 
vations or a space station (ISchutz \911\ . In addition, there is 
the renewed interest of major space agencies for Lagrangian 


point colonization. Eurthermore, de Fillipi 1978 has made a 
review of the ideas of \0’Neill & Gerard 7974l about build¬ 
ing space colonies at the /L 4 and L 5 positions. These space 
stations could be used as a waypoint for travel to and from 
the cislunar space. Currently, there are no probes orbiting 
L 4 or L 5 points for any celestial couple. |5a/azar et al. 2012 \ 
proposed an alternative transfer strategy to send spacecraft 
to stable orbits around the Lagrangian equilibrium points L 4 
and L 5 based on trajectories derived from the periodic orbits 
around L 5 . In the present work or main goal is to provide a 
method for keeping a station around the orbit of L 5 . 

In Section 2 we describe our models and intial condi¬ 
tions. In Section 3 we present the results of the simulations 
as phase portraits in the vicinity of L 5 and give the escape 
rate and the associated average lifetime of the test parti¬ 
cles in the different models. The orbit-keeping procedure 
is demonstrated with the use of a test case and described 
in detail in Section 4. Einally our conclusions are given in 
Section 5. 


2 Model and initial conditions 


In the present study two models were used; the (re¬ 


stricted) three body-problem (R3BP) (Szebehely 19671 and 


(restricted) four body-problems (R4BP). In our simulations 
one of the bodies, the spacecraft has negligible mass, but 
we have kept its mass in the equations of motion in order 
to be able to compute the fuel consumption (hence the re¬ 
stricted adjective is in parentheses). The R3BP models the 
motion of a particle under the gravitational attraction of two 
larger bodies (also called primaries), where the primaries 
are point masses that revolve in elliptic orbits around their 
common centre of mass. In our case the primaries were the 
Earth and the Moon. In order to take into account of the ef¬ 
fect of the Sun the model was extended to the R4BP. The 
equations of motion are valid in 3D and were transformed 
to the geocentric reference frame. The model is depicted in 


Z 



Figure 1 The model in the inertial reference frame 
OXYZ: R, Re, Rm, Rs denotes the position vectors of 
the spacecraft. Earth, Moon and the Sun respectively. The 
vectors r, and r 2 denote the geocentric positions of the 
spacecraft. Moon and the Sun respectively. 


Fig.m where R, Re, Rm, Rs denotes the position vec¬ 
tors of the spacecraft. Earth, Moon and the Sun in an inertial 
reference frame. The relative vectors r = (x, y, z) ri and 
r 2 denote the geocentric position vector of the spacecraft. 
Moon and the Sun respectively. 

The first three Lagrangian equilibrium points, Li,L 2 
and L 3 are aligned with the two larger bodies. On the other 
hand, L 4 and L 5 points form equilateral triangles with the 
two masses on the plane of orbit whose common base is the 
line between their center of mass and L 4 leads the orbit of 
the smaller mass and Li follows (Eig.|2]). 

Of particular interest is the stability of the Lagrangian 
points; the three collinear stationary points are always un¬ 
stable, whereas the stability of L 4 and L 5 points depends on 
the mass ratio between the two larger bodies ( |Danby 1964 1 
and eccentricity. In the case of the important celestial cou¬ 
ples like Earth-Moon or Sun-Earth L 4 and L 5 points are sta¬ 
ble. This property makes the fuel required for a spacecraft 
to maintain its relative position there to be almost zero. 

The first goal was to map the phase space around the L 5 
point in the RTPB of the Earth-Moon system. To explore 
the stable and unstable regions of the motion of the space¬ 
craft in the R3BP model, the equations of motion have to be 
integrated for a large set of initial conditions. Eor the given 
initial condition of the Moon (see Table[T]i, we integrated or¬ 
bits of the test particle placing it at different points as start¬ 
ing points along the line connecting L 5 and the Earth as 
demonstrated in Eig.|2l The initial velocity was chosen such 
that the angular velocity of the test particle was equal to that 
of the L 5 point. The time of integration was 1300 days (for 
several test cases it was extended to 144 000 days). The de¬ 
termination of the length of the integration time is explained 
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Figure 2 Initial conditions: top view of the geocentric 
system. The initial positions were chosen from between the 
light gray bars located on the line connecting the Earth and 
the Ls point. The notion are the same as in Fig. [T] 


in Section 3.1. This initial condition generation procedure 
was applied also in the simulations of the R4BP (see Fig. 

0 . 

The major goal of the present study is to present a proce¬ 
dure to keep a probe on orbit around the EMLsQ point of the 
Earth-Moon system. The initial condition for the satellite is 
given in Table [T] 

We applied a variable step-size Runge-Kutta method 
(IDormand & Prince 1978l l. The initial conditions of the 
major bodies are from the JPL database with epoch = 
2455999.5, in geocentric ecliptic Cartesian coordinate sys¬ 
tem, see Table [U 


3 Results 

The elliptic restricted three-body problem is often used in 
dynamical investigations of planetary systems or to de¬ 
sign transfer strategy to send spacecrafts to stable orbits 
around one of the Lagrangian points. However, in real mis¬ 
sions other perturbations must be taken into account. In the 
present study we include the most relevant force, which 
is the Sun’s gravitational attraction. Fig. [3] clearly demon¬ 
strates the crucial role the Sun plays in the orbital evolution 
of a probe: in the R3BP the probe moves on a regular orbit, 
whereas in the R4BP its trajectory is chaotic and the space¬ 
craft escapes in less than 1 year. 

* The five Sun-Earth Lagrangian points are called SEL1-SEL5, and 
similarly those of the Earth-Moon system EML1-EML5, etc. 
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Figure 3 The xi (t) of the test particle in the R3BP (con¬ 
tinuous line) and in the R4BP (dashed line). The test particle 
was started from the vicinity of the L 5 point of the planar 
R3BP Earth-Moon system. 


First we mapped a larger region of the phase space 
around L 5 in different models. The aim was to acquire a 
general picture of the structure of the phase space: where 
are the stable and chaotic regions and what are their sizes 
and shapes. Also we were interested in the escape rates and 
in the average lifetime of the test particles. This knowledge 
is essential to design the orbit and the orbit-keeping strategy 
of a probe. 


3.1 Transient chaos and escape rates in the vicinity of 

Lb 


In transient chaos, typical trajectories, i.e. trajecto¬ 
ries initiated from random initial conditions, escapes 
any neighborhood o f the non attracting chaotic set 
dKovacs & Erdi 20091 . To define the escape rate one dis¬ 
tributes a large number of initial points Nq at time fp in 
a phase space region R and then follows the trajectories of 
these points. A quantity measuring how quickly particles 
leave R is the escape rate k\ 


N{t) = 


( 1 ) 


where N{t) is the number of particles which are still part of 
the initial region R. A small value of k implies weak ’’re¬ 
pulsion” of typical trajectories by the non attracting chaotic 
set. The escape rate is a property solely of the non attracting 
chaotic set. 

To find K for the R3BP using the initial conditions de¬ 
scribed in the previous section the space region R was de¬ 
termined. Using a regular grid covering R the trajectories 
of the test particles were followed for 1 000 days both in 
the R3BP and R4BP. The number of test particles we used 
for the R3BP is 400 and for the R4BP is 10 000. Our main 
model was the R4BP, we investigated it first and we used 
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Model 

A 

B 

T[day ] 

R3BP 

5.9858 

-0.00224152 

446.126 

R4BP 

7.7848 

-0.00324845 

307.839 


Table 2 Results of the linear fit to the models R3BP and 
R4BR In the first column the model in the second the con¬ 
stant A and the third column displays the slope of the fitted 
linear model. The escape times are given in the last column. 


more initial points. After receiving the result, (Fig. 4b), we 
saw that for the R3BP will be enough less initial points to 
get the escape rate, because it is defined only by the slope 
of the line. The results are shown in Fig. |4]left and right 
panel, respectively. The average lifetime r is estimated as 
T « 1/k. We have fitted paired data to the 

linear model, y = A + Bx, by minimizing the chi-square 
etTor statistic. The slope of the line is the parameter k and 
its reciprocal is the average lifetime. The average lifetime 
of the test particles are approximately 440 days (ft! 16 rev¬ 
olution of the Moon) and 310 days (ft 11 revolution of the 
Moon) in R3BP and R4BP, respectively. The later is sig¬ 
nificantly shorter than that in R3BP. The parameters of the 
fitted line and t are summarized in Table |2] 


According to the results obtained for the average life¬ 
time the simulation time was set as usual the 4-6 xr 
dTel & Gruiz 20061 ILai & Tel 201 II) . i.e. in all our simula¬ 
tions the test particles were integrated for T = 1300(ft 5r) 
days. The results are shown in Fig. |5] In the figures the 
horizontal axis is the test particle’s initial x coordinate in 
units of the average distance tq = 400000 km of the 
Earth-Moon. The vertical axis is the test particles initial x 
velocity in rg/day unit. The number of the test particles 
was N = 10® in the case of Fig. |5] panel a and b, and 
N = 4 X 10® for panel c and d. 

In the first row of Fig. |5] the (x, x) € 
[-1.014,-0.7142] X [-0.2128,0.21082] region was 
mapped. First the planar case of the R3BP was studied 
(Fig. |5]panel a). It is clearly visible that - as expected - the 
region around L 5 is stable (black zone) and there exits other 
stable regions as black area. The motion started exactly 
from L 5 is regular as it is well visible from Fig. [3 The x{t) 
of a test particle started exactly from the L 5 point (see the 
large cross in Fig.|5]panel a) is shown in Fig. [3 


The size of the quasi-periodic region next to point L 5 
without the Sun (panel a) coincides with the stable re gion 
determined by a different method by \Erdi et al. 2009\ The 
largest distance of the initial points of the stable region from 
L 5 in the outward direction of the Earth is around 0.005 rg 


which is almost the same found bv lErdi et al. 2009 


In the next set of runs we have included the Sun and 
studied the orbits in the planar case. The most striking dif¬ 
ference is the almost total disappearance of the stable region 
around L 5 . The perturbation from the Sun destroys most of 
this stable zone and renders this region to a strongly chaotic 
area of the phase space. A small island remains on the left 
of the L 5 . It is worth noting that other larger stable regions 


survived the perturbing effects of the Sun and remained sta¬ 
ble with approximately the same structure and extension. A 
dashed line bounding box shows the region which was stud¬ 
ied with more initial particles, that region is shown in panel 
c. 

In the second row of Eig.|3we have explored the closer 
vicinity of L 5 with larger resolution (we have kept the 
ranges of the axes fixed to ease the comparison of the re¬ 
sults). The reason to zoom in the phase space was to have a 
much finer resolution of the stable and chaotic regions in the 
vicinity of L 5 . We have also increased the number of initial 
conditions to 4 million (the model for panel c was the same 
as that of panel b). In panel c a smaller plot shows a stable 
region shifted to the left of L 5 . This region is stable for the 
integration time, but will vanish for longer times. 

In the last panel of Eig.|3the computation was extended 
to the spatial case. One can see an extremely complex struc¬ 
ture and also it useful to note that the escape rate is de¬ 
creased with the introduction of the 3rd dimension. 


4 Chaos control 


Control of chaos refers to a process wherein a tiny pertur¬ 
bation is applied to a chaotic system, in order to obtain 
a desirable (chaotic, periodic, or stationary) behavior. The 
idea of chaos control was published at the beginning of this 
decade by Ott, Grebogi & Yorke 1990 the ideas for control¬ 
ling chaos were outlined and a method for stabilizing an 
unstable periodic orbit was suggested, as a proof of prin¬ 
ciple. The basic idea consisted in waiting for a passage of 
the chaotic orbit close to the desired periodic behavior, and 
then applying a small fittingly chosen perturbation, in order 
to render an the chaotic motion more stable and predictable, 
which is often an advantage. 


According to the method developed by 
\Ott, Grebogi & Yorke 7990| (hereafter OGY method), 
first one obtains information about the chaotic system by 
analyzing a Poincare section of the phase space. After 
the information about the section has been gathered, the 
system is integrated until it comes near a desired periodic 
orbit in the section. Next, the system is forced to remain 
on the target orbit by perturbing the appropriate parameter. 
One strength of this method is that it does not require 
a detailed model of the chaotic system but only some 
information about the Poincare section. It is for this reason 
that the method has been so successful in controlling a wide 
variety of chaotic systems (e.g. turbulent fluids, oscillating 
chemical reactions, electronic circuits). 


There exits other alternative methods for chaos control. 
In general the strategies for the control of chaos can be clas¬ 
sified into two main classes, namely; closed loop or feed¬ 
back methods and open loop or non feedback methods. 

The first class includes those methods which se¬ 
lect the perturbation based upon a knowledge of the 
state of the system. Among them, we mention (besides 
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Figure 4 The escape rates for the models R3BP (left panel) and R4BP (right panel). The continuous curve is the measured 
data, the dashed line is the fitted line N{t) = A + Bt. The parameters of the line A and B as well as t = —XjB are 
displayed on each panel (see also Table|2|i. 


\Ott, Grebogi & Yorke 1990) the so called occasional pro¬ 
portional feedback simultaneously introduced hv \Hunt 1991\ 
and \Petrov et al 1993\ and the method introduced by 
\Pyragas 1992\ All these methods are model independent, 
in the sense that the knowledge on the system necessary to 
select the perturbation can be done by simply observing the 
system for a suitable learning time. 

The second class includes those strategies which con¬ 
sider the effect of external perturbations (independent on the 
knowledge of the actual dynamical state) on the evolution 
of the system. Periodic (ILima & Pettini 1990t or stochastic 
dBraiman & Goldhirsch 199T1 ) perturbations have been seen 
to produce drastic changes in the dynamics of chaotic sys¬ 
tems, leading eventually to the stabilization of some peri¬ 
odic behavior. These approaches, however, are in general 
limited by the fact that their action is not goal oriented, i.e. 
the final periodic state cannot be decided by the operator. 

4.1 Orbit-keeping procedure 

In our case the parameters of the system consist only the 
mass of the celestial bodies, therefore it is not possible to di¬ 
rectly apply the OGY method to control the chaos. A differ¬ 
ent approach will be discussed and we stress that the method 
is based on observation and numerical experiments. 

Chaos is a double-edge sword, since it destabilizes the 
motion and renders it unpredictable, but at the same time it 
offers the opportunity to transport spacecrafts to distant lo¬ 
cations with minimal fuel consumption. The present method 
is the result of an extensive numerical experiment. The idea 
was that if we cannot apply the OGY method, i.e. ’’kick” the 
body onto the stable manifold of the hyperbolic point, then 


let us try to drive it away from that point. In order to accom¬ 
plish this a small force (thrust) was added to the equation of 
motion of the test particle: 

Fo iff G [fci, tci + Ati], . ^ 

0 otherwise, ’ * ~ ' 

where Fq > 0 is the magnitude of the applied force and 
is the unit vector in the direction of the particle’s velocity, 
tci is the time instant when the control begins and Ati is the 
length of the control. 

As a proof of our concept we will show a test case where 
we have applied our orbit-keeping procedure for a space¬ 
craft to keep it around the L 5 point. In Fig.| 6 ]we have plot¬ 
ted the time evolution of the x coordinate of the spacecraft 
(its initial conditions are in the last column of Table [T]|. In 
the upper panel the probe orbits without control and stays 
around L 5 for approximately 300 days, then it escapes that 
region and begins to oscillate widely. In the lower panel the 
x{t) is plotted when the above defined force (thrust) was ap¬ 
plied from tci = 66.32 days for Ati = 29.87 days (see the 
vertical lines), which corresponds to about one period of the 
Moon. The value of Fq = 9.01 x 10“® was used in this case. 
This value is the result of dozens of tests and proved to be 
a quite general value and not only the ’’best” for the present 
situation. We have also investigated braking force {Fq < 0) 
but this alternative in our tests never stabilized the orbit. 

To determine the ith start time of the control td one has 
to integrate the motion of the probe and monitor its orbit. 
The x{t) function and the x—y and x—x sections are proved 
to be useful in order to determine the start and stop time of 
the control. Our numerical experience has shown us that the 
control must start well before the escape and at a moment 
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Figures The x — x phase space. The points corresponds to orbits which did not leave the x C [—2, 2] and x € [—2, 2] 
region after 1300 days. The legend shows the dimension of the simulation (2D or 3D), whether the Sun was included or 
not, (R3BP or R4BP) and the number of initial points (JV). The position of the L 5 point is shown by a large cross. Panel a) 
shows the result of the planar R3BP; b) shows the result of the planar R4BP; c) shows a smaller region around L 5 in the 
planar R4BP, while d) shows the result of the spatial R4BP for 4 million test particles. 


when the orbit intersects the cc — y plane, i.e. when the z = 0 
and i > 0 conditions are met. To find the ’’best” starting 
time is the result of a systematic search but the structure or 
pattern of the points on the x — y and x — x section provides 
some additional information. In Fig. |2]the x — y section with 
z = 0 and i > 0 of the test particles trajectory from f = 0 to 
t < 3000 days is shown. The orbit without control is repre¬ 
sented by red diamonds while the one with control is shown 
by blue plus signs. The first few section points are numbered 


to see the time evolution of them. In Fig.0the x — x section 
with y = 0 and y > 0 of the test particles trajectory from 
f = 0 to f < 3000 days is demonstrated. The blue contin¬ 
uous line segments connects the successive section points 
of the trajectory with control, the red dotted line segments 
connects the points without control. The controlled trajec¬ 
tory’s section points stay inside a bounded region whereas 
those without control leave this region after a few period of 
the Moon. 
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Figure 6 The x{t) of the probe in the R4BP without (up¬ 
per panel) and with control (lower panel). The probe was 
started from the vicinity of the L 5 point (see Table [T]l. The 
orbit without control stays around L 5 for approximately 300 
days, then it leaves that region. We have applied a thrust 
from tci = 66.32 days for Afi = 29.87 days. This control 
stabilizes the motion up to approximately 5000 days. On the 
lower panel the vertical lines denote the start and stop of the 
control. 


In the present example shortly after the first section 
point the sequence of points do not follow any ’’regular” 
pattern, (i.e. a ’’scalloped edge” ellipse or a snowflake on 
the X — y section plane, see the blue plus signs in Fig.|2]or 
a distorted ellipse on the x — x section plane, see the blue 
continuous line in Fig. [ 8 ]l therefore we switch on the addi¬ 
tional force immediately at point No. 3. There are two other 
parameters which we have to set: (/) the force Fq and {ii) 
the length of the control . 

We have experiment with the magnitude of the force Fq 
and found that the value of Fq w 9.01 x 10“® could be 
applied in all the cases we have tried. 

In all our cases the start and stop time of the control 
is a time instant of a section in the x — y section plane. 
In our test cases the stop time was either the next, the sec¬ 
ond or the third section after the start time, so the duration 
of the control Ati is equal to the time span between two or 
more subsequent points on the x — y section plane. The time 
span between two successive section point is approximately 
equal to the period of the Moon, therefore loosely speaking 
we may measure the length of the control in ’’periods”. One 
period of a control means the time difference between two 
successive section points. We emphasize that these results 
based solely on numerical observations and experimenta¬ 
tion: we have extensively studied other Af^s but it turned 





- 1.0 - 0.5 0.0 0.5 1.0 

X [rj 

Figure? The x — y section with z = 0 and i > 0 of 
the test particles orbit from f = 0 to f < 3000 days. The 
red diamond symbols are the section points of the trajec¬ 
tory without control. The blue plus sign symbols are the 
section points of the trajectory with control. The numbers 
over the symbols are the nth section point. No 1, 2 and 3 
points are overlapping, since the control starts at point No. 
3. The Earth is in the centre. 

out that the best choice was always one of the moment of a 
successive section. 

The control in the test case was started at the 3rd point 
on the x — y section plane (see Fig. |7]i which corresponds 
to time tci = 66.32 days and the thrust was applied until 
the next section was detected, i.e. for one period. Having 
switched off the thrust we have integrated the system and 
monitored the orbit of the spacecraft. As it can be seen form 
Fig.|9]the x{t) function oscillates quasi periodically about 0 
and its amplitude does not varies significantly {x{t) in Fig. 
|9]is a continuation of x{t) in Fig. | 6 ]l. This behavior can be 
observed till t ~ 5000 days when the probe leaves the vicin¬ 
ity of L 5 again. Just as in the first case, we have applied a 
thrust from fc 2 = 4364.28 days until we have detected the 
next section, which resulted in At 2 = 26.6 days (somewhat 
shorter than the in previous control). The same value of Fq 
was used again. 

The trajectory on the x — y and x — x section planes 
are shown in Fig. [TO] and Fig. (TT] respectively. In Fig. [TO] 
we see the continuation of the pattern that was observed in 
Fig. 12] and in Fig. [TT] we see the continuation of Fig. [ 8 ] In 
Fig. [To] the x — y section with z = 0 and i > 0 of the 
test particles trajectory from t — 3000 days to f < 6000 
days is shown, the symbols are the same as in Fig. [7] The 
red diamonds and the blue plus signs overlap until when 
the control was switched on and then the points depart from 
each other. The start time and length of the control is again 
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Figure 8 The x — x section with y = 0 and y > 0 of the 
test particles orbit from f = 0 to f < 3000 days. The red 
dotted line connects the section points of trajectory without 
control, while the blue continuous line connects those with 
control. The numbers are the nth section point of the trajec¬ 
tory without control. No 1 and 2 are overlapping, since the 
control starts at point No. 3. 
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Figure 10 The x — y section with z = 0 and i > 0 of the 
test particles orbit from t = 3000 days to f < 6000 days. 
The symbols are the same as in Fig. [T] The Earth is in the 
centre. 

a result of a systematic search: tc2 and At 2 were found to 
result in a quasi periodic orbit for f = 18 900 days. In Fig.fTTI 
the X — X section with y = 0 and y > 0 of the test particles 
trajectory from t = 3000 days to f < 6000 days is shown. 
The controlled trajectory’s section points (blue continuous 
line) stay inside a bounded region whereas those without 
control (red dotted line) leave this region after a few period 
of the Moon. We note that the size of the bounded region 
hardly changed comparing it with the previous time span 
shown in Fig. [ 8 ] The result of the systematic search was that 
the start time of the second control was tc 2 = 4364.28 days 
and the length of the control was one period (At 2 = 26.6 
days). 

We used a third and a fourth control beginning at t 
= 18899.49 days and t= 18979.37 days which resulted in 
bounded motion until t = 30196.72 days. We note that the 
duration of the third control was again one period but three 
periods for the fourth control. The results in terms of the 
x{t), X — y and x — x are very similar to the ones presented 
above, so we do not show them. 

The orbit-keeping procedure can be summarized as fol¬ 
lows: 


Figure 9 The x{t) of the test particle. The orbit after the 
first control without the next control stays around L 5 for 
approximately 5000 days, then it leaves that region. In the 
lower panel we have applied a thrust from tc 2 = 4364.28 
days for At 2 = 26.6 days. This control stabilizes the motion 
up to approximately 18 900 days. The thick vertical lines 
denote the start and stop of the control. 


1. Determine the moment of ’’escape”: the motion of the 
probe is integrated until it leaves the region allowed by 
some prescription (i.e. the mission instructions). 

2. A systematic search is carried out in order to find the 
best or acceptable time instant tci when the control 
should start. According to our numerical observations 
(0 this time instant should correspond to the moment 
of a section in the x — y plane and iii) it must precede 
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Figure 11 The x — x section with y = 0 and y > 0 of the 
test particles orbit from t = 3000 days to t < 6000 days. 
The symbols are the same as in Fig. 

with a few to a few dozens of section points the time of 
escape. 

3. Another systematic search is carried out in order to hnd 
the length of the control Atf. according to our obser¬ 
vation the stop of the control is also determined by the 
moment of section on the x — y plane; in our cases this 
section point was the first, the second or the third after 
the start time. 

4. We have found that Fq ~ 9 x 10“® was applicable in all 
our test cases. 

This process was tried for several initial conditions and 
we could stabilize the motion of the spacecraft in all test 
cases: the probe was temporary pushed on a bounded or¬ 
bit. In the demonstrated case the time intervals of the quasi 
periodic motion were long enough for a usual space mis¬ 
sion to carry out. After the hrst control it was approximately 
13 years, then 38 years and after the last control 31 years. 
We have repeated our test cases with probes having differ¬ 
ent masses: the start time of the control and its duration 
(29.87H-26.6 days that resulted quasi periodic motion for 13 
+ 38 years) was the same for probes as massive as 10® kg. 
This is understandable since the most massive probe’s in¬ 
fluence on the dynamics of the other bodies of the model is 
still negligible. We have estimated the energy needed for the 
control: for a one period length control of 1 kg mass satellite 
the consumed energy is approximately 12 kJ. 

5 Conclusion 

We have applied the concept of SMD and the basic idea 
of the OGY method to stabilize the motion of a spacecraft 


around the L 5 in the R4BP. The phase space structure in 
the vicinity of L 5 was mapped in both the R3PB and in the 
R4BP. We have studied the escape rate of the trajectories 
and fitting a line to the measured data we were able to esti¬ 
mate the average lifetime of the particles which turned out 
to be 446 (R3BP) and 307 (R4BP) days. Integrating the two 
different models we mapped the phase space around the L 5 
point of the Earth-Moon system. 

Using the phase space portrait our main goal was to ap¬ 
ply the OGY method to keep a spacecraft close to the vicin¬ 
ity of L 5 . The OGY method cannot be directly applied to 
our problem at hand (there are no variable system param¬ 
eters) therefore we modified the equation of motions with 
the addition of a time dependent force to the motion of the 
spacecraft. There are three free parameters of our orbit¬ 
keeping procedure: (/) the magnitude of the thrust (Fq), (ii) 
the start time of the control and {in) the length of the con¬ 
trol. We have performed extensive numerical experiments 
and according to our observations we found the following 
rule of thumb useful: (/) Fq ~ 9 x 10“® is applicable in 
all cases, (ii) the start time of the control is a time instant 
of a section in the x — y section plane and {Hi) the length 
of the control is an integer multiply of the period (the time 
between two successive section points). 

We note that we lack the physical background and 
explanation of our orbit-keeping procedure. We are fully 
aware of the fact that the underlying physics and the proper¬ 
ties of the chaos responsible for the observed behavior must 
be carefully studied in the future. We are performing sys¬ 
tematic research in this direction in order to understand the 
basic processes causing the observed behavior. 
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Moon 

Sun 

Spacecraft 

X 

-0.5166166629896271 

368.66644026568872 

-0.9141820107443692 

y 

-0.7573376905382426 

-47.60083686893512 

0.0687343088982397 

z 

-0.0172345994799509 

0.0004280643569276 

0 

X 

0.185382724767311 

0.9306857451293241 

-0.0252733218674244 

y 

-0.1362138844144604 

6.4039891664344024 

-0.2286530912785001 

z 

0.0197096442582322 

0.0000507123227653 

0 


Table 1 The initial condition of the bodies taken from the JPL. The values are valid for JD = 2455999.5. The initial 
condition for the spacecraft is in the close vicinity of L 5 . 
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